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Abstract — We  consider  the  problem  of  cognitive  code-division 
channelization  (simultaneous  power  and  code-channel  allocation) 
for  secondary  transmission  links  co-existing  with  an  unknown 
primary  code-division  multiple-access  (CDMA)  system.  We  first 
develop  a  blind  primary-user  identification  scheme  to  detect  the 
binary  code  sequences  (signatures)  utilized  by  primary  users.  To 
create  a  secondary  link  we  propose  two  alternative  procedures 
-one  of  moderate  and  one  of  low  computational  complexity-  that 
optimize  the  secondary  transmitting  power  and  binary  code¬ 
channel  assignment  in  accordance  with  the  detected  primary 
code  channels  to  avoid  “harmful”  interference.  At  the  same 
time,  the  optimization  procedures  guarantee  that  the  signal- 
to-interference-plus-noise  ratio  (SINR)  at  the  output  of  the 
maximum  SINR  linear  secondary  receiver  is  no  less  than  a  certain 
threshold  to  meet  secondary  transmission  quality  of  service  (QoS) 
requirements.  The  extension  of  the  channelization  problem  to 
multiple  secondary  links  is  also  investigated.  Simulation  studies 
presented  herein  illustrate  the  theoretical  developments. 

Index  Terms — Blind  user  identification,  code-channel  allo¬ 
cation,  code-division  multiple-access,  cognitive  radio,  dynamic 
spectrum  access,  power  allocation,  signal-to-interference-plus- 
noise  ratio. 


I.  Introduction 

WITH  the  rapid  proliferation  of  a  variety  of  consumer 
oriented  wireless  devices,  demand  for  access  to  radio 
spectrum  has  been  growing  dramatically  and  the  limited 
available  spectrum  is  becoming  increasingly  congested.  At 
the  same  time,  location-dependent  bands  of  pre-licensed  radio 
spectrum  may  experience  low  utilization  [1],  Cognitive  radio 
(CR)  is  an  emerging  technology  aiming  at  improving  spectrum 
utilization  efficiency  by  allowing  secondary  users/networks  to 
opportunistically  share  radio  spectrum  originally  licensed  by 
primary  users/networks  without  causing  “harmful”  interfer¬ 
ence  to  them  [2]-|5]. 

Cognitive  radio  networks  can  be  categorized  according  to 
two  modes  of  operation:  cooperation  mode  and  coexistence 
mode  [6],  In  cooperation  mode,  primary  users  cooperate 
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with  secondary  users  and  share  information  to  avoid  mutual 
interference.  In  coexistence  mode,  there  is  no  form  of  coop¬ 
eration  and  secondary  users  must  have  the  ability  to  detect 
the  presence  of  primary  users  [7],  [8]  and  change  behavior 
accordingly  to  avoid  mutual  interference. 

Past  work  in  the  young  field  of  cognitive  code-division 
channelization  includes  coexistence  power  control  1 9]  as  well 
as  distributed  resource  allocation  of  spectral  bands,  power, 
and  data  rates  among  multiple  secondary  users  for  multi¬ 
carrier  CDMA  systems  [10].  Cooperation-mode  bit  rate  and 
spreading  factor  adjustments  for  a  secondary  CDMA  system 
under  interference-minimizing  code  assignments  were  carried 
out  in  [11].  In  [12],  [13],  a  secondary  maximum  signal- 
to-interference-plus-noise  ratio  (SINR)  code-division  link  is 
designed  subject  to  SINR  requirements  for  the  primary  system 
which  is  presumed  known  (cooperation-mode  cognitive  radio). 
Outside  the  framework  of  cognitive  code-division  altogether, 
interesting  work  in  the  form  of  joint  beamforming  and  power 
allocation  algorithms  was  reported  in  [14]-[18].  In  particular, 
in  [14]  the  radio  frequency  spectrum  of  interest  was  divided 
into  a  set  of  multiple  orthogonal  channels  and  was  shared 
between  primary  and  secondary  networks  using  orthogonal 
frequency  division  multiple  access  (OFDMA).  In  [15]-[17], 
joint  spatial-channel  and  power  allocation  algorithms  for  cog¬ 
nitive  radio  networks  were  developed.  In  [18],  the  authors 
provide  a  solution  for  leasing  spectrum  for  a  fraction  of  time 
to  secondary  users  based  on  the  idea  that  secondary  nodes  can 
earn  spectrum  access  in  exchange  for  transmission  assistance 
to  the  primary  link  (cooperative  communications  paradigm). 

In  this  work,  we  focus  on  coexistence-mode  cognitive  radio 
and  investigate  the  problem  of  establishing  a  secondary  code¬ 
division  link  coexisting  with  an  unknown  primary  CDMA 
system.  In  particular,  we  investigate  -to  the  best  of  our 
knowledge  for  the  first  time  in  the  context  of  cognitive  radio- 
the  problem  of  blindly  identifying  the  binary  codes/signatures 
utilized  by  primary  users  when  neither  channel  state  infor¬ 
mation  nor  pilot  signaling  (training  sequence)  is  available. 
Then,  we  study  the  design  of  a  power  and  binary-code- 
channel  allocation  protocol  for  the  secondary  link  that  will 
not  cause  “harmful”  interference  to  the  existing  primary  users. 
Since  post-processing  interference  sensing  is  not  feasible  in 
coexistence  mode  cognitive  radio,  to  quantify  “harm”  we 
use  the  periodic-total-square-correlation  (PTSC)  interference 
metric  in  our  optimization  problem  as  the  mathematical  means 
to  protect  co-channel  primary  users  [2],  [19],  [20].  At  the 
same  time,  to  satisfy  quality-of-service  (QoS)  requirements  for 
the  secondary  link,  the  power  and  code-division  optimization 
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problem  is  constrained  to  have  S1NR  at  the  output  of  the 
maximum  SINR  linear  receiver  of  the  secondary  link  no  less 
than  a  certain  threshold.  We  recognize  that  the  above  described 
fundamental  cognitive  code-division  radio  formulation  is,  re¬ 
gretfully,  a  non-convex  NP-hard  problem.  Yet,  using  herein 
existing  SINR-maximization  signature  design  methodologies 
we  are  able  to  develop  novel,  realizable  suboptimum  so¬ 
lutions  of  varying  computational  complexity  with  excellent 
cognitive  system  performance  characteristics  as  demonstrated 
by  simulation  studies  included  in  this  paper.  The  theoretical 
developments  and  experiments  can  be  readily  extended  to 
cover  multiple  secondary  links  alongside  the  primary  CDMA 
system. 

The  rest  of  this  paper  is  organized  as  follows.  The  code¬ 
division  cognitive  radio  problem  formulation  is  presented  in 
Section  II.  A  novel  primary  user  identification  algorithm  is 
introduced  in  Section  III.  The  power  allocation  and  signature 
design  procedure  for  a  secondary  link  is  developed  in  Section 
IV.  In  Section  V,  we  extend  the  procedure  to  solve  for  multiple 
secondary  links.  Simulation  results  are  presented  in  Section  VI 
and,  finally,  a  few  conclusions  are  drawn  in  Section  VII. 

II.  System  Model  and  Problem  Formulation 

The  following  notation  is  used  throughout  this  paper.  Bold¬ 
face  lower-case  letters  indicate  column  vectors  and  boldface 
upper-case  letters  indicate  matrices;  C  denotes  the  set  of  all 
complex  numbers,  ()T  and  ( )H  denote  the  transpose  and 
transpose-conjugate  operation,  respectively;  is  the  L  x  L 
identity  matrix,  DTe{-}  denotes  the  real  part  of  a  complex 
number,  sgn{-}  denotes  zero-threshold  quantization,  and  E{-} 
represents  statistical  expectation.  Finally,  |  •  |,  ||  •  ||,  and  ||  •  ||F 
are  the  scalar  magnitude,  vector  norm,  and  matrix  Frobenius 
norm,  respectively. 

In  the  following,  we  consider  a  primary  code-division  sys¬ 
tem  with  a  primary  transmitter  PT  and  K  primary  receivers 
PRi,i  =  1 , 2, ....  A,  as  shown  in  Fig.  1.  The  primary  trans¬ 
mitter  (for  example,  base  station)  PT  communicates  downlink 
with  the  K  primary  receivers  PRi,  i  1,2,...,  A,  over  dis¬ 
tinct  code-division  channels  defined  by  individual  normalized 
binary  codes/signatures  s =  -^={±l}i,  i  =  1,2,..., A, 
where  L  is  the  signature  length  (system  processing  gain).  We 
consider  also  a  potential  concurrent  secondary  code-division 
link  in  the  spectrum  band  of  the  primary  downlink  channel 
between  a  secondary  transmitter  ST  and  receiver  SR.  If 
the  primary  system  is  frequency-division-duplex  (FDD),  then 
the  secondary  ST-to-SR  link  will  work  on  the  frequency 
licensed  to  the  downlink  channel  of  the  primary  system  and 
can  operate  at  any  time.  If  the  primary  system  is  time-division- 
duplex  (TDD),  then  the  secondary  link  can  transmit  data  only 
during  the  downlink  slot  of  the  primary  system.  The  secondary 
communication  link  is  activated,  whenever  possible,  with  a 
(normalized)  binary  signature  c  =  and  transmitting 

power  P  >  0.  All  transmitted  signals,  primary  and  secondary 
when  appropriate,  are  assumed/modeled  to  propagate  over 
multipath  Rayleigh  fading  channels  and  experience  additive 
white  Gaussian  noise  (AWGN). 

We  first  assume  that  the  secondary  transmitter  ST  is  quiet 
and  we  examine  how  the  signal  sent  by  PT  is  observed  by 
SR.  After  carrier  demodulation,  chip  matched  filtering  and 


Fig.  1.  Primary/secondary  code-division  system  model  of  a  primary  trans¬ 
mitter  PT,  K  primary  receivers  PRi,i  =  1,2 and  a  secondary 
transmitter/receiver  pair  ST,  SR  (all  received  signals  exhibit  multipath 
Rayleigh  fading). 


sampling  at  the  chip  rate  over  a  presumed  multipath  extended 
data  bit  period  of  LM  =  L  +  M  —  1  chips  where  M  is 
the  number  of  resolvable  multipaths,  the  observed  data  vector 
y (n)  £  CLm  by  SR  takes  the  following  general  form1 

K 

y{n)  =  Y,  y/Pjbi(n)  Hsj  +  i  +  n,  n  =  l,2, ...,  (1) 

i= 1 


where  H  £  CtMXi  is  the  multipath  channel  matrix  between 


PT  and  SR 

hi 
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.  0 
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.  0 
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.  0 
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with  entries  hrn 

£  C, 

TO  =  1,  . 

..,M, 

considered  as 

com- 

plex  Gaussian  random  variables  to  model  fading  phenomena, 
\fEl  >  0  and  6,(n)  £  {±1}  are  the  amplitude  level  and  nth 
transmitted  bit  of  primary  user  i,  i  =  1 , ,K,  respectively, 
i  £  CLm  denotes  multipath  induced  inter-symbol-interference 
(ISI),  and  n  is  a  zero-mean  additive  white  Gaussian  noise 
(AWGN)  vector  with  autocorrelation  matrix  <r2l£M.  The 
information  bits  6,;(n)  are  viewed  as  binary  equiprobable 
random  variables  that  are  independent  within  a  user  stream 
(in  n  =  1,2,.. .)  and  across  users  (in  z  =  1,2,...,  K).  Since 
the  effect  of  ISI  is,  arguably,  negligible  for  most  applications 
of  practical  interest  where  the  number  of  resolvable  multipaths 
is  much  less  than  the  processing  gain,  for  mathematical  con¬ 
venience  we  will  not  consider  the  ISI  terms  in  our  theoretical 
developments  that  follow2.  Thus,  the  primary  users’  signal 

Equation  (1)  assumes  that  SR  is  symbol-synchronous  to  PT.  While  this  is 
not  a  technical  requirement  at  all  for  the  secondary  link  design  problem  set  and 
solved  herein  based  on  the  periodic  cross-channel  total  squared  correlation, 
it  greatly  simplifies  the  presentation  and  notation  of  this  material. 

-However,  naturally  ISI  will  be  considered  and  accounted  for  in  our 
simulation  studies. 
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observed  by  SR  in  (1)  is  simplified/approximated  by 

K 

y (n)  =y^v/^'&4(n)Hs4  +  n,  n  =  1,2,....  (3) 

i- 1 

In  our  cognitive  system  model,  the  secondary  link  is  taken 
to  be  chip-synchronous  to  the  primary  network  (worst  case 
interference  scenario)  with  the  same  chip  rate  and  symbol- 
synchronous  (see  Footnote  1).  Also,  without  loss  of  generality 
and  for  simplicity  in  notation,  we  assume  that  the  multipath 
channels  between  PT  and  PRi,  i  =  1 ,K,  PT  and  SR, 
and  ST  and  SR,  all  have  the  same  number  of  resolvable  paths. 
Then,  when  the  secondary  communication  link  is  activated 
with  a  (normalized)  binary  signature  code  c  =  -^={±1}L  and 
transmit  power  P  >  0,  the  aggregate  signal  vector  received 
by  SR  can  be  expressed  as 

r(n)  =  VPb(n)Gc  +  y(n),  n=  1,2,...,  (4) 

where  G  €  CLmxL  is  the  ST  to  SR  channel  matrix  with 
multipath  channel  coefficients  grn  G  C,  m  =  1, . . . ,  M,  and 
y(n)  is  given  by  (3). 

Information  bit  detection  at  SR  is  carried  out  via  linear 
maximum  SINR  filtering  (or,  equivalently,  minimum  mean 
square  error  filtering)  as  follows 

b{n)  =  sgn {53e{w,^aa;S7Wflr(n)}}  ,  n  =  1,2,...,  (5) 

where  w maxSiNR  =  cR-1Gc  G  CLm,  c  >  0,  is 

the  maximum  SINR  filter  and  R  =  E{yyff}  is  the  au¬ 
tocorrelation  matrix  of  the  y(n)  signal  in  (3)  that  con¬ 
stitutes  primary-system  disturbance  to  SR.  Practically,  R 
is  estimated  by  averaging  over  N  >  LM  observation 
samples  r (n)  when  ST  is  silent  (P  =  0),  R(iV)  = 
7jltn=ir(n)r(ri)H  =  ±  £^=i  y(n)y  {n)H .  The  output 
SINR  of  the  filter  w maxSiNR  can  be  calculated  to  be 

T  4  E{lWmaxS7jVfl(vAP^Gc)|2}  =  1G(,  (fi) 

®  {  \^maxSINRy\^  } 

To  attain  a  certain  QoS  level  for  the  secondary  link,  we  need 
to  jointly  design  the  binary  signature  c  and  the  transmitting 
power  P  to  have  the  SINR  value  T  at  the  output  of  the  max¬ 
imum  SINR  filter  w maxsiNR  no  less  than  a  given  threshold 
7  >  0,  i.e.  T  >  7. 

At  the  same  time,  due  to  the  coexistence  of  the  secondary 
link  with  the  primary  network,  interference  is  introduced  to 
the  primary  receivers.  The  secondary  link  is  to  be  allowed  to 
activate  only  if  the  interference  to  each  primary  receiver  is  not 
“harmful.”  The  difficulty  is  that  in  cognitive  radio  networks 
operating  in  coexistence-mode,  the  primary  network  does  not 
cooperate/talk  to  secondary  users  and  the  latter  do  not  have 
global  knowledge  of  system  parameters,  such  as  the  multipath 
channel  coefficients  [hi 7, . . . ,  between  PT  and  PRi,  or 
the  channel  coefficients  [<7,^, . . . ,  'q^m]  between  ST  and  PRi, 
i  =  1,  2, . . . ,  K  (see  Fig.  1),  or  the  primary  binary  signatures 

i  =  1, . . . ,  K,  in  (3)  and  the  filters  utilized  by  each  primary 
user  receiver3.  Therefore,  post-processing  interference  sensing 
is  not  feasible. 

3  In  fact,  even  if  global  system  parameter  values  (code  channels,  transmis- 
sion  power  values,  channel  coefficients,  receive  filters)  were  to  be  provided 
to  the  secondary  system  by  the  primary  system,  SINR  optimized  secondary 
link  design  under  primary-system  SINR  constraints  is  a  non-convex,  NP-hard 
problem  [12],  [13]. 


Toward  a  realistic,  realizable  solution  to  the  problem  of 
coexistence  secondary  link  design,  we  recall  that  the  periodic 
(cyclic)  total  squared  cross-correlation  (PTSC)  value  [19],  [20] 
is  a  useful  measure  to  evaluate  multiple  access  interference 
(MAI)  when  channels  exhibit  multipath  behavior.  In  this  spirit, 
we  propose  to  use  the  PTSC  value  as  a  metric  to  evaluate 
the  interference  caused  by  the  secondary  link.  For  notational 
simplicity,  let  s,|;,  i  =  1  denote  the  cyclic  right- 

shifted  version  of  s4  G  -^{±1}L,  by  l  bit  positions  when 
l  =  0, 1, . . . ,  L  —  1,  and  cyclic  left-shifted  version  of  s by  l 
bit  positions  when  l  =  0,  —1, . . . ,  1  —  L,  (hence,  s.j|0  =  s,). 
The  PTSC  between  signatures  c  and  S  j  for  multipath  shifts  up 
to  lag  M  is  defined  as 

M 

PTSC(c, Sj)  =  ^  \cTSill\2,i=l,...,K.  (7) 

l=—M 

In  this  context,  we  define  the  generalized  correlation  in¬ 
terference  by  a  secondary  link  with  power  P  >  0  to  the  ith 
primary  receiver  as 

Ii  =  P  ■  PTSC(c,  sQ,  i  =  1, . . .  ,K.  (8) 

We  understand  that  'L,  serves  as  a  simple  “worst-case”  measure 
of  the  effect  of  the  secondary  link  on  the  ith  primary  user. 
Then,  in  this  context,  the  secondary  link  can  be  activated  by 
assigning  a  signature  c  and  power  P  >  0  if  the  interference  to 
every  primary  user  is  less  than  a  threshold  Ith  >0:1,;  <  Zth 
Vz  =  1  If 

Tmax  =  max{li  :  i  =  1, . . . ,  K}  (9) 

is  the  strongest  generalized  interference  to  primary  receivers, 
the  activation  condition  is  equivalent  to  Zmax  <  Zth- 

Our  objective  is  to  jointly  design  the  binary  signature  c 
and  the  transmitting  power  P  >  0  for  the  secondary  link  to 
minimize  Imax  under  the  constraint  that  the  secondary  link 
achieves  its  pre-determined  SINR  requirement  7: 

minimize  Tmax  —  maa;{P-PTSC(c,  s,)  :  i  =  1, . . . ,  K} 
p> 0,  ce^{±i}L 

(10) 

s.  t.  r  =  PctGhR_1Gc  >  7.  (11) 

Then,  if  the  resulting  minimized  Imax  is  indeed  less  than  Tth, 
the  secondary  link  can  be  activated;  otherwise,  it  is  kept  idle. 

The  cognitive  code-division  channelization  problem  formu¬ 
lated  in  (10),  (11)  requires  knowledge  of  all  active  primary- 
user  binary  signatures  s ,,  i  =  l .....  K ,  for  the  evaluation 
of  the  created  interference  in  the  form  of  PTSC(c,  s,), 
i  =  1, ...  ,K.  However,  in  our  assumed  coexistence  mode 
of  operation,  the  secondary  network  cannot  obtain  such 
knowledge  directly  from  primary  networks.  Therefore,  before 
starting  to  solve  the  channelization  problem  in  (10),  (11),  the 
secondary  network  needs  to  blindly  detect  (i)  the  number  of 
active  primary  users  K  and  (ii)  their  binary  signatures  s(, 
i  =  1 .....  A .  With  respect  to  item  (i)  (population  size  iden¬ 
tification  problem),  we  can  utilize  for  example  the  algorithm 
developed  recently  in  [21].  Due  to  space  limitations,  in  this 
paper  we  do  not  deal  further  with  this  issue  and  instead  assume 
that  K  is  correctly  identified.  With  respect  to  item  (ii),  in  the 
next  section  we  develop  an  iterative-least-square  (ILS)-based 
procedure  that  can  blindly  detect  the  primary  users’  binary 
signatures  s,,  i  =  1, . . . ,  K,  from  the  observed  primary  signal 

y  (n). 
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TABLE  I 

Iterative  least-squares  procedure 


1)  d  :=  0;  initialize  B^0)  E  {dzl}iCxAr  arbitrarily. 

2)  d  :=  d  +  1; 


\Hd)  :=  Y(B(d))T 

g(d)  .=  sgn|<Hc  [(v(d_1) )H (V(d_1) )j  1  (v(d-1))ffY  j. 

3)  Repeat  Step  2  until  (BW.VW)  =  (bG-F.vG-1)). _ 


III.  Primary-system  Identification 

If  we  denote  the  (energy  inclusive)  channel  processed 
signature  by 

v^y/Elnsi,  i  =  l,2,...,K,  (12) 

then  the  observed  signal  in  (3)  can  be  expressed  as 

K 

y(n)  =  Vjbi{n)  +  n  =  Vb(n)  +  n,  n  =  l,2, ...,  (13) 

i- 1 

where  V lm*k  —  [vi , . . . ,  v/c]  is  the  effective  signature 
matrix  and  b(n)  =  [&i(n), . . . ,  6/f(n)]T  is  the  vector  of  bits 
for  all  K  users  at  the  nth  transmission  period.  If  SR  is  able 
to  collect  N  observation  vectors  y (n),  n  =  1, 2, . . . ,  N,  then 
(13)  can  be  rewritten  in  matrix  form  as 

Y  =  VB  +  N  (14) 

where  Y  £  CLmxW  is  the  observation  matrix,  B  = 
[b(l), . . . ,  b(iV)]  is  the  K  x  N  data  matrix  that  contains  the 
N  bits  transmitted  for  each  of  the  K  primary  users,  and  N  is 
an  Lm  x  N  Gaussian  noise  matrix. 

To  detect  the  binary  signatures  s,;,  i  =  1, ....  K,  we  first 
estimate  the  channel  processed  signature  set  V  from  the 
observation  matrix  Y.  Our  approach  begins  by  formulating 
the  signature  set  estimation  problem  as  a  joint  detection  and 
estimation  problem  with  the  following  least  squares  (LS) 
solution 

V,  B  =  arg  min  ||Y-VB||*.  (15) 

Be{±i}(**V  "F 

x  k 

The  above  LS  solution  is  maximum-likelihood  optimal  as 
long  as  N  is  white  Gaussian.  In  any  case,  regretfully,  joint 
detection  and  estimation  by  (15)  has  complexity  exponential 
in  N K .  We  consider  this  cost  unacceptable  and  attempt  to 
reach  a  quality  approximation  of  the  solution  by  alternating 
least  squares  estimates  of  V  and  B,  iteratively,  as  described 
below. 

A.  Iterative  Least  Squares  Procedure 

The  basic  idea  behind  such  an  iterative  least  squares  (ILS) 
solution  [22] -[25]  is  to  compute  an  LS  update  of  one  of 
the  unknown  (matrix)  parameters  conditioned  on  a  previously 
obtained  estimate  of  the  other  (matrix)  parameter  and  continue 
on  until  convergence  is  observed.  The  iterative  least-squares 
procedure  for  the  solution  of  our  problem  in  (15)  is  presented 
in  Table  I.  Superscripts  in  Table  I  denote  the  iteration  index. 
Derivation  details  are  presented  in  the  Appendix. 


We  understand  that  convergence  of  the  developed  iterative 
least  squares  procedure  to  a  globally  optimal  LS-solution  of 
(15)  is  not  guaranteed  in  general.  The  quality  (least-squares 
fit)  of  the  end  convergence  point  depends  heavily  on  the 
initialization  point  and  arbitrary  initialization  -which  at  first 
sight  is  unavoidable  for  blind  primary  system  identification- 
offers  little  assurance  that  the  iterative  scheme  will  lead 
us  to  appropriate,  “reliable”  (close  to  minimal  least-squares 
fit)  solutions.  Re-initialization  and  re-execution4  of  the  ILS 
procedure  is  always  possible  but  the  challenge  is  how  to  assess 
whether  solutions  returned  by  the  ILS  procedure  are  reliable 
or  not  without  any  side  information  or  pilot  signaling.  The 
rest  of  this  section  is  devoted  to  addressing  this  challenge. 

Since  B  and  V  are  jointly  detected  and  estimated,  corre¬ 
spondingly,  if  one  is  not  reliable,  neither  is  the  other  in  general. 
We  first  examine  the  reliability  of  the  bit  matrix  decision 
B  =  |  b  | . . . . ,  b^-]7  returned  by  the  ILS  procedure  of  Table  I. 
The  sample  cross-correlation  between  any  two  bit  streams  is 

Vi,j  =  bf hj/N,  i^j,i,j  =  l,...,K.  (16) 

Formally,  the  true  information  bits  are  independent  within 
user  streams  and  across  users.  If  j  were  to  be  viewed 
as  approximately  normally  distributed  with  zero  mean  and 
variance  then  the  probability  of  | %,j|,  i  ^  j,  being 
larger  than,  say,  the  threshold  value  -)=  is  very  low  at 
about  0.3%  (we  can  calculate  >  ^=)  «  0.003). 

Motivated  by  this  calculation,  we  introduce  below  Criterion 
1  that  classifies  convergence  points  of  the  ILS  procedure  in 
Table  I  as  “unreliable”  based  on  the  sample  statistics  of  the 
returned  data  matrix  B. 

Criterion  1:  If  \r]ij  \  >  for  some  i  ^  j  £  (1,2,...,  K}, 

then  (B,  V)  returned  by  the  ILS  procedure  in  Table  I  are 
classified  as  “ unreliable .”  ■ 

Criterion  1  provides  the  means  for  coarse  identification  of 
unreliable  solutions.  An  unreliable  convergence  point  would 
then  trigger  re-initialization  and  re-execution  of  the  ILS 
procedure  in  Table  I.  To  enhance  the  end  accuracy  of  the 
blind  primary-system  identification  procedure,  we  propose  one 
additional  criterion  based  on  the  returned  estimated  effective 
signature  matrix  V.  We  will  motivate  our  proposal  by  exam¬ 
ining  the  normalized  cross-correlation  between  the  estimated 
channel  processed  signatures  v,  returned  by  the  ILS  procedure 
upon  convergence  and  the  true  channel  processed  signatures 
*  =  1, ... ,  K.  Based  on  N  =  128  snapshots  for  a  system 
with  K  =  4  primary  users  with  equal  6dB  SNR  and  processing 
gain  L  =  31  (Gold  signature  codes),  we  run  the  Criterion 
1-equipped  ILS  procedure  20  times.  The  distribution  of  the 
twenty  returned  solutions  of  0,  =  p^-|  ^  y  for  signal  i  =  1 
in  Fig.  2  (representative  of  all  other  signals)  reveals  that  (i) 
Criterion  1  is  not  by  itself  sufficient  to  eliminate  erroneous 
solutions,  however  (ii)  there  exist  “reliable”  regions/clusters 
in  which  most  of  the  Criterion  1-equipped  ILS  convergence 
points  lie  close  to  the  true  channel  processed  signatures.  The 
basic  idea  behind  our  second  and  final  refinement  of  the  ILS 
blind  primary-system  identification  procedure  is  to  identify 
and  average  these  reliable  clustered  estimates.  Of  course, 

4In  practical  implementations,  re-initialization  and  re-execution  of  the  ILS 
procedure  may  also  be  needed  whenever  numerical  instabilities  create  a  rank 
deficient  matrix  BO)  in  Table  I  and  bG)bG)t  is  not  invertible. 
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TABLE  II 

Cross-Correlation  Enhanced  ILS 


For  j  :=  1  to  D 

1)  Execute  ILS  of  Table  I  with  arbitrary  initialization  and  obtain 
estimates  v,:.  i  =  1, . . . ,  K. 

2)  If  estimates  are  reliable  according  to  Criterion  1, 

let  :=  v^,  i  =  1, .  . . ,  K: 
else  go  to  1). 

End 

For  i  1  to  K 

3)  Identify  reliable  estimates  of  v;  according  to  Criterion  2. 

4)  Calculate  average  over  all  reliable  estimates  v,  by  (19). 

End 


5)  Set  V  =  [vi, . . . ,  vjf], 

6)  Execute  ILS  of  Table  I  with  initialization 
/— if— \  -1  —H 


§(°)  =  sgn  91c 


V  V 


V  Y 


Fig.  2.  Experimentation  with  ILS  of  Table  I  under  Criterion  1 :  Distribution 
of  normalized  cross-correlation  between  vi  and  vi  after  twenty  runs  (K  = 
4,  L  =  31,  N  =  128,  SNR  *  =  6dB,  i  =  1,2, 3, 4). 

identification  of  the  reliable  estimates  is  not  a  trivial  task  due 
to  our  complete  lack  of  knowledge  of  v,;  (or  s,). 


initialized  at  B(°)  =  sgn  <  Die 


rv 


V  Y 


We  dub 

ILS  with  both  Criteria  1  and  2  Incorporated  Cross-Correlation 
Enhanced  ILS  (CC-ILS)  and  summarize  the  complete  proce¬ 
dure  in  Table  II. 


In  this  context,  assume  that  we  have  D  estimates  of  v,; 
denoted  by  vp'P  i  =  1 , ...  ,K,  j  =  1 ,D,  obtained  by 
D  runs  of  the  Criterion  1 -equipped  ILS  procedure.  From  the 
example  of  Fig.  2,  we  notice  that  reliable  estimates  vp)  of 
Vi  have  high  normalized  cross-correlation  (close  to  1)  with 
each  other,  while  they  have  low  normalized  cross-correlation 
with  other  unreliable  estimates  of  Vj.  In  contrast,  unreliable 
estimates  tend  to  have  low  normalized  cross-correlation  with 
each  other.  Therefore,  the  reliability  of  vp*  may  be  quantified 
by  examining  the  sum-cross-correlation  with  the  rest  vpP  t  ^ 


j. 


D 

=  E 

t=M/j 


K  Vi  I 
llvpplllvfpl  ' 


(17) 


A  reasonable  threshold  value  for  binary  reliability  classifica¬ 
tion  may  be  the  average  value 


. ^  d8) 

j= 1 

utilized  in  the  proposed  Criterion  2  below. 

Criterion  2:  Let  vp'  be  the  estimates  of  v,  returned  by 
D  arbitrary  initializations  of  the  Criterion  1 -equipped  ILS 
procedure  of  Table  I,  i  =  1, . . . ,  K,  j  =  1, . . . ,  D.  If  >  ~pi, 
then  vp^  is  considered  a  reliable  estimate  of  the  v^;  otherwise 
we  declare  it  as  unreliable.  ■ 

Next,  we  average  our  reliable  (according  to  Criterion  2) 
estimates  of  the  channel  processed  signatures  v,  to  produce 
one  last  high-quality  initialization  of  the  ILS  algorithm  of 
Table  I.  Let  Si  denote  the  set  of  all  reliable  estimates  of  v,; 
according  to  Criterion  2  and  let  |tS,  |  denote  the  cardinality  of 
Si.  Our  averaged  estimate  of  the  matrix  V  is  now  given  by 
V  with 


V  = 


where  v, 


1 

N 


E 

j&Si 


i  =  l...,K, 
(19) 


i.e.  Vi  is  the  average  over  all  reliable  estimates  of  v,  according 
to  Criterion  2.  We  execute  ILS  in  Table  I  a  final  time 


B.  Detection  of  Primary  Binary  Signatures 

After  obtaining  an  estimated  (energy  inclusive)  channel 
processed  signature  set  V  by  CC-ILS  of  Table  II,  we  develop 
another  procedure  to  extract  the  individual  primary  binary 
signatures  s i  =  by  decomposition  of  V.  The 

channel  processed  signatures  can  be  rewritten  as 


vi  =  v^Sih,i  =  l,...,A-,  (20) 

where  h  =  [hi, ... ,  Hm]T  and 

f  s4l)  0  1 


Si  = 


Si(L) 


Si(l) 


(21) 


0 


Si(L) 


LmxM 


If  the  binary  signatures  s  i,  i  =  1, ...  ,K,  were  known,  by  (20) 
we  could  estimate  \/Tfh  as 

y/Eih  =  (Sf Sj)-1Sp v^  i  =  l,...,K,  (22) 

where  v,  is  the  ith  column  of  matrix  V.  Then,  a  quality 
estimate  of  h  (scaled)  could  be  produced  by  averaging, 

~  1  K  - 

h  =  ^E^  (23) 

i=l 

to  create  a  matrix  channel  estimate  H  by  (2).  Given  H  (and 
Vi),  one  could  detect  the  binary  signatures  by 

s,  =  -^sgn{me{(HHH)-1HHVi}},i  =  l,...,A7  (24) 

The  proposed  individual  binary  signature  extraction  algorithm 
from  V  is  now  ready.  Initialize  s i  e  -pj={±l}L,  i  =  1, . . . ,  K , 
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Fig.  3.  Performance  of  blind  primary-users  identification  algorithm:  Primary 
transmitter  PT  has  K  downlink  users  with  length  L  —  31  Gold  signatures. 
The  transmitted  signal  propagates  over  a  3 -path  Rayleigh  fading  channel  to 
secondary  receiver  SR  that  collects  N  observation  samples  to  run  Table  II 
procedure  followed  by  execution  of  (22)-(24). 


arbitrarily,  and  alternate  computation  between  (22),  (23),  and 
(24),  iteratively.  Stop  when  convergence  is  observed5. 

To  illustrate  briefly  the  proposed  primary  binary  signature 
detection  algorithm,  we  consider  a  primary  downlink  CDMA 
system  with  K  =  4  or  I\  =8  active  users  that  utilize 
Gold  signatures  with  length  L  =  31.  The  primary  users’ 
signals  are  transmitted  with  equal  power  over  a  multipath 
Rayleigh  fading  channel  with  M  =  3  resolvable  paths  in 
the  presence  of  additive  white  Gaussian  noise.  SR  is  able 
to  collect  N  =  256, 384,  or  512  observation  samples  and 
employs  our  proposed  method  to  extract  the  primary  users’ 
binary  signatures.  The  experiment  is  repeated  105  times  with 
randomly  drawn  channel  coefficients.  In  Fig.  3,  we  plot  the 
probability  of  correct  identification  of  all  binary  signatures, 
Pr(si  =  SjVi  =  1,.. . , K),  as  a  function  of  SNR.  It  can  be 
seen  that  even  under  low/moderate  SNR  values  the  proposed 
method  can  correctly  extract  all  signatures  with  sufficient 
sample  support  (secondary  receiver  observation  time  interval). 
It  is  worth  pointing  out  that,  experimentally,  when  errors  do 
occur  only  one  or  two  signatures  have  few  chip-bit  errors  only. 

Having  addressed  the  blind  primary  system  identification 
problem  -arguably-  satisfactorily,  in  the  next  section  we  at¬ 
tempt  to  solve  the  secondary  power  and  code-channel  alloca¬ 
tion  problem  of  (10)  to  enable  most  frequent  activation  of  the 
secondary  link. 


IV.  Secondary  Code-Division  Link  Channelization 


The  secondary  link  has  to  satisfy  the  SINR  constraint  (11); 
namely,  that  the  SINR  of  the  secondary  receiver  SR  at  the 
maximum-SINR  linear  filter  output  is  no  less  than  a  QoS 
requirement  7:  P crGffRGc  >  7,  7  >  0.  To  that  respect, 
it  suffices  to  set  the  power  value  P  >  0  to 


P  = 


7 

cTGffRGc 


(25) 


sThe  iterative  procedure  may  be  re-executed  with  distinct  initialization  if 
convergence  cannot  be  observed  after  sufficient  iterations. 


and  the  SINR  constraint  (11)  is  always  satisfied  with  equality. 
Then,  the  optimization  problem  in  (10)  is  equivalent  to 


minimize  Tmax  =  max{P  ■  PTSC(c,  s ,)  :  i  =  1, . . . ,  K } 

s.  t.  P=  ctg2RGc-  (26) 

The  maximum  generalized  correlation  interference  Imax  is 
the  product  of  two  components,  transmit  power  P  and  maxi¬ 
mum  PTSC  value.  Thus,  to  minimize  Imax  we  need  to  design 
a  binary  signature  c  to  minimize  the  product  of  the  required 
transmit  power  P  =  ctgi?RGc  times  max  {PTSC(c,  s^)}. 

This  problem  is  still  non-convex  NP-hard  (see  Footnote  1  and 
[12],  [13]).  We  must,  therefore,  pursue  (disjoint)  suboptimal 
design  procedures  if  we  wish  to  keep  the  computational 
complexity  manageable.  At  first,  we  look  at  minimizing  P 
alone.  The  binary  signature  c  that  minimizes  P  maximizes 
the  denominator  of  (25)  which  is  the  SR  output  SINR  with 
unit  ST  transmit  power: 


c  =  arg  max  c1  Ac 


(27) 


where  A  =  G^RG.  At  this  point,  the  SINR-maximizing 
binary  signature  designs  of  polynomial  complexity  developed 
in  [26],  [27]  can  be  used  directly.  We  recall  that  in  [26] 
the  binary  signature  vector  is  optimized  under  a  rank-2  ap¬ 
proximation  of  the  matrix  A,  while  in  [27]  the  arcs  of  least 
SINR  decrease  from  the  real  maximum  SINR  solution  are 
evaluated.  Both  algorithms  first  generate  L  candidate  binary 
signatures6,  denoted  by  qj  £  -4^{±1}L,  j  =  1, . . . ,  L,  which 
can  provide  high  output  SINR.  Then,  the  signature  among 
them  with  highest  output  SINR  is  selected.  While  the  highest- 
SINR  signature  minimizes  the  transmitting  power  P  required 
to  satisfy  any  given  QoS  constraint,  it  may  result  to  high  values 
of  PTSC  with  respect  to  primary  user  signatures  and  conse¬ 
quently  let  the  secondary  link  introduce  strong  interference  to 
the  primary  system.  Therefore,  we  propose  to  evaluate  all  L 
binary  signatures  qj  £  ^={±1}L,  j  =  1, ...  ,L,  returned  by 
the  solver  of  (27)  in  [26]  or  [27]  in  our  interference  metric 
for  the  v'th  primary  user 


=  Pj  ■  PTSC(q,, Si),  j  =  1, . . . , L-  i  =  1, . . . ,  K,  (28) 

and  find  the  maximal  generalized  correlation  interference 
value  caused  by  qj,  Ijtmax  =  maxjlj^  :  i  =  1,  ...,K}, 
j  =  1,...,L.  Then,  we  choose  the  signature,  power  pair 
(qj» ,  Pj » )  which  has  the  least  maximal  interference.  If  the 
resulting  maximal  interference  by  (qr  ,  Pj* )  is  introduced 
to  the  i*th  primary  user,  i.e.  Xj*,i*  =  min{2j  max  :  j  = 
1, . . . ,  L},  and  still  Zj.^.  <  Zth,  the  secondary  link  is  allowed 
to  access  the  channel  by  assigning  signature  c  =  qy  -  and 
power  P  =  Pj . .  We  refer  to  this  method  of  selecting  a  pair 
of  signature  and  power  from  the  candidate  set  as  passive 
interference  suppression  and  outline  the  procedure  in  Table 
III. 


The  secondary  link  design  method  in  Table  III  exclusively 
focuses  on  transmit  power  minimization.  The  PTSC  factor 

6  We  recall  that  [27]  finds  ( LT  —  T  f  1 )  binary  sequences  that  are  closest  to 
T  arcs  of  least  SINR  decrease  in  the  I2  sense.  In  this  paper,  we  only  consider 
one  slowest  descent  arc  and  generate  L  binary  sequences.  This  is  sufficient 
to  closely  approximate  the  performance  level  reached  when  all  L  —  1  slowest 
descent  arcs  are  considered. 
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TABLE  III 

Secondary  Link  Signature  and  Power  Design  w /  Passive 
Interference  Suppression 


TABLE  IV 

Secondary  Link  Signature  and  Power  Design  w/  Active 
Interference  Avoidance 


Input  A  :=  G"RG 

Obtain  q j,  j  =  1, . . . ,  L  as  solution  candidates  for  (27)  by  [26]  or  [27]. 
Calculate  Pj,  j  =  1, . . . ,  L,  by  (25). 

Calculate  Xj,i,  j  =  1, . . . ,  L,  i  =  1, . . . ,  K,  by  (28)  and 

Xj,max  —  TflCLX \Xj } i  •  %  —  1,  .  .  .  ,  K}  • 

Select  j*,  i*  such  that  Xj*  ,i*  =  mm{Zj?max  :  j  =  1, . . . ,  L}. 

Output  P,  i* ,  qj* ,  and  Xj  *  ,* . 

<Ith, 

transmit  on  channel  qj*  with  power  Pj*', 
else  seize 


is,  of  course,  evaluated  and  accounted  for  in  (28)  but  not 
actively  optimized  (minimized).  To  further  reduce  the  max¬ 
imum  interference  caused  to  the  primary  system  and  further 
improve  the  chances  of  spectrum  sharing,  we  next  propose 
to  iteratively  adjust  the  binary  channel  signature  to  actively 
avoid  interference  by  jointly  reducing  the  PTSC  value  with 
the  most  impacted  at  each  time  i*th  primary  user’s  signature 
,  as  well  as  minimize  the  transmit  power. 

Define  S,:.  =  [s;.|_Mr  ■  L,  s*. |0,  •  •  • ,  s^m]  and  calcu¬ 
late  PTSC(c,  s,. )  =  cTSi*Sf,c.  To  combine  the  SINR- 
optimization  problem  of  (27)  with  the  PTSC  suppression 
task,  after  executing  the  procedure  in  Table  III  we  update 
A  =  G;/RG  —  aS,*  S'C  where  a  >  0  is  an  introduced 
weighting  factor.  Then,  by  re-executing  the  procedure  in 
Table  III  with  the  updated  A,  we  obtain  a  new  optimal 
pair  (q  j*,Pj*)  and  new  maximum-interference  Ij*,i* .  If  the 
new  maximum-interference  1j*,i*  is  reduced,  we  iteratively 
update  A  4—  A  —  o  S,  -  Sf*  and  re-execute  the  algorithm  in 
Table  III.  This  procedure  will  be  stopped  when  the  maximum 
interference  cannot  be  suppressed  further.  We  outline  the 
procedure  in  Table  IV. 

Duplexing  of  the  secondary  link  can  be  implemented  in 
TDD  mode  by  switching  the  roles  of  transmitter  and  receiver. 
The  reverse  secondary  link  can  be  activated  by  a  new  pair  of 
transmit  power  and  channel  code/signature  values  calculated 
by  the  same  algorithm  in  Table  III  or  IV  to  satisfy  reverse 
link  QoS  requirements.  The  following  section  is  devoted  to 
the  generalization  of  the  channelization  problem  to  multiple 
secondary  receivers  (one-to-many  secondary  downlink  trans¬ 
missions). 

V.  Design  of  Multiple  Secondary  Code-Division 
Links 


Input  R,  G,  and  s;,  L<=  1, . . . ,  K 
d  :=  0. 

A  :=  GflRG. 

Obtain  Pj*,  q j*,  and  Xj*,i*  by  Table  III. 

Pd  :=  Pj* ,  cd  :=(ij*,I^ax  '.=  Ij*ti*. 

While  d  =  0  or  Xi ax  <  X^ 
d  :=  d  +  1; 

A  <—  A  —  aSi-Sf,; 

obtain  Pj*,  q j*,  and  Xt*  ;*  by  Table  III; 

pd  p  rd  —  jd  j 

r  ,—  j-j*,^  . —  ,  -i,max  ■ —  -Lj*  jt*  . 

End 

Output  P(d_1) ,  c(d_1),  and  Imax^- 
If  Imo  <  2Jth' 

transmit  on  channel  cV'-  1 '  with  power  p(d~ O; 
else  seize. 


Y 


ra. 


PT 


ST 


PPk 


SR, 


SR„ 


Fig.  4.  Primary/secondary  code-division  system  model  of  a  primary 
transmitter  PT,  K  primary  receivers  PRi,i  =  1,2,...,  if,  a  secondary 
transmitter  ST,  and  Q  secondary  receivers  SRi,  i  =  1,  2, . . . ,  Q  (all  receives 
signals  exhibit  multipath  Rayleigh  fading). 


The  discrete  time  received  signal  vector  of  the  fcth  SR  can 
be  expressed  as 


In  this  section,  we  extend  the  study  to  cover  code-division 
channelization  for  multiple  secondary  links.  Particularly,  as 
shown  in  Fig.  4,  we  consider  the  scenario  where  an  ST 
attempts  to  communicate  downlink  with  Q  potential  SRs7 .  Let 
C  =  {1, . . . ,  Q}  be  the  set  of  all  secondary  link  indices  which 
attempt  to  share  the  radio  spectrum.  Among  the  Q  potential 
secondary  links,  Ks  >  0  secondary  links  will  be  activated  by 
assigning  to  each  a  (normalized)  signature  Ck  =  -^{±1}L 
and  transmitting  power  Pf~  >  0,  e5C  C  where  S  represents 
the  set  of  active  secondary  links. 

7  Multiple  one-to-many  secondary  downlinks  may  be  treated  sequentially 
establishing  one  secondary  downlink  at  a  time  using  techniques  presented  in 
this  section,  treating  previously  established  secondary  downlinks  as  part  of  a 
"virtual  primary  system." 


Yk{n)  =  sfl\bk{n)GkCk+  ^  i/P^bj(n)GkCj  +  yk(n), 
k£S,n  =  1,2,...,  (29) 

where  Gfc  £  CLMXi  is  the  channel  matrix  constructed  by 
the  multipath  channel  coefficients  \gk)M,  ■  ■  ■ ,  9k, m]  between 
the  secondary  transmitter  and  the  fcth  secondary  receiver  and 
yk  represents  comprehensively  MAI  from  all  primary  users 
plus  AWGN.  The  autocorrelation  matrix  of  yk  is  denoted  by 
Rfc  =  Ejyfcyj^}.  The  output  SINR  of  the  fcth  secondary 
receiver  filter  -wmaxSiNR,k  is  given  by 

Tk=Pk  c^GfR^GfeCfc  (30) 
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where 


Rfc  =  Gfc  ^  Pjcjcj  Gf  +  Rfc  (31) 

\jes,j^k  ) 

is  the  autocorrelation  matrix  of  the  combined  channel  dis¬ 
turbance.  To  satisfy  a  certain  QoS  requirement,  each  active 
secondary  link  is  supposed  to  surpass  a  pre-determined  SINR 
value  7^  >  0,  i.e.  we  need  Tk  >  "fk  for  each  k  £  S.  At  the 
same  time,  interference  caused  to  the  primary  system  needs  to 
be  kept  in  check.  To  evaluate  the  interference  introduced  by 
all  active  secondary  links  to  the  zth  primary  user,  we  define 
the  cumulative  generalized  correlation  interference 

I,:  =  5ZPfe-PTSG(Cfc’Si)’  *  =  1>- (32) 

fees 

Our  objective  is  to  maximize  the  number  of  active  sec¬ 
ondary  links  Ks  =  |«S|  under  the  constraints  that  each 
secondary  link  achieves  its  pre-determined  SINR  7 k  and  the 
cumulative  interference  to  each  primary  user  is  less  than  a 
threshold  2th  >  0: 


maximize  Ks 

s.t.rfe(cfc,Pfc)  >7 Vfc  €  S  and 2,;  <2th  Vi  =  1, . . . ,  K.(33) 


This  problem  is,  again,  non-convex  NP-hard.  Motivated  by  the 
single  secondary  link  channelization  algorithm  of  Table  IV 
(with  active  interference  avoidance)  in  the  previous  section, 
we  propose  instead  to  -centrally  at  the  secondary  transmitter 
(base  station)  ST-  successively  activate  one  secondary  link  at 
a  time  in  a  way  that  it  introduces  conditionally  least  cumulative 
interference  until  the  accumulated  interference  exceeds  the 
threshold  Xtll  for  some  primary  user.  The  details  of  the 
algorithm  are  described  as  follows.  Initially  set  S  =  0. 
Utilizing  the  algorithm  in  Table  IV,  we  individually  eval¬ 
uate  each  secondary  link  k  £  C  and  obtain  the  signature 
c k,  transmitting  power  l\,  and  the  maximum  generalized 
interference  value  2 max,k  caused  to  a  primary  user  by  the 
potential  secondary  link  with  design  (c k,Pk)-  Among  the  Q 
potential  links,  we  select  the  one,  say  the  /th  secondary  link, 
which  introduces  least  maximum  generalized  interference, 
i.e.  2 max,t  —  rnin{Imax,k  •  k  £  C}.  If  Tmax,t  <c  2th, 
the  fth  secondary  link  is  turned  on  with  code  channel  ct 
and  power  Pt\  if  link  t  did  go  through,  we  continue  on 
examining  the  remaining  C  —  {t }  candidates  by  the  algorithm 
in  Table  IV  under  the  cumulative  generalized  interference 
metric  in  (32)  since  one  secondary  link  is  already  activated. 
Say  potential  secondary  link  2  £  C  —  {i}  is  the  one  with 
max,z  =  tnin\Tmax,k  •  k  £  C  {V}}  and  still  2maX)2  ^  Pth 
(if  Tmax,z  f_  2th,  link  z  is  rejected  and  no  further  study  is 
required).  Since  the  incoming  secondary  link  z  also  interferes 
with  the  existing  secondary  link  t  and  may  degrade  its 
output  SINR  below  the  minimum  acceptable  level,  we  adjust 
iteratively  the  transmitting  power  of  each  secondary  link  by 
calculating 


Pk 


Tk 


cTGPH71Gkck 


k  =  t,z, 


(34) 


until  convergence  is  observed.  With  the  updated  powers,  we 
evaluate  the  cumulative  generalized  interference  in  (32)  one 
more  time.  If  all  2 ,,  z  =  1, . . . ,  K ,  are  below  the  threshold,  the 
zth  secondary  link  is  turned  on  and  we  move  on  to  identify 


TABLE  V 

Cognitive  Code-Division  Channelization  for  Multiple 
Secondary  Links 


C  :=  (1, . . .  ,Q};  S  :=  0. 

While  |C|  >  0 

Calculate  R*.  VfcgC  by  (31). 

Obtain  Pk,  ck  TmaXtk  Vfc  £  C  by  Table  IV. 

Find  t.  e  C  s.  t.  lmax,t  =  min{Imax  k  :  k  £  C}. 

Update  powers  Pk,  k  £  (5,  t}  by  (34). 

Calculate  Zj,  i  =  1, . . . ,  K,  by  (32)  (with  the  updated  powers). 
If  Z;  <Zth  V-i  =  1,. ..  ,K, 

C  =  C\{t}-  S  =  S\J{t}. 
else 

C  =  0. 

Endif 

End 

Output  {( Pj ,  Cj)  :  j  £  5}  (activated  secondary  channels). 


additional  secondary  links  with  the  same  procedure;  otherwise, 
z  is  rejected  and  no  more  secondary  link  are  allowed  to  share 
the  spectrum.  The  complete  algorithm  is  summarized  in  Table 
V. 


VI.  Experimental  Simulation  Studies 

To  illustrate  the  presented  algorithmic  developments,  we 
consider  a  primary  multiuser  CDMA  system  with  I\  active 
users  that  utilize  Gold  signatures  with  length  L  =  31.  The 
primary  users’  signals  are  transmitted  with  equal  per-user 
power  Ei  =  E2  =  ...  =  Ek  —  lOdB  over  a  multipath 
Rayleigh  fading  channel  with  M  =  3  resolvable  paths  in 
the  presence  of  additive  white  Gaussian  noise.  At  first,  one 
secondary  link  attempts  to  share  the  spectrum  with  target 
receiver  filter  output  SINR  7  =  lOdB. 

The  secondary  link  code  channel  and  power  are  optimized 
by  the  algorithm  in  Table  III  (passive  interference  suppression) 
and  Table  IV  (active  interference  avoidance).  Specifically, 
we  utilize  the  rank-2  SINR-maximization  binary  signature 
design  method  [26]  and  set  the  weighting  factor  for  the 
algorithm  in  Table  IV  to  a  =  0.05.  If  the  resulting  maximum 
generalized  correlation  interference  is  less  than  the  threshold 
2th,  the  secondary  link  is  allowed  to  share  the  spectrum.  The 
simulation  experiment  is  repeated  105  times.  With  varying 
primary  user  population  K  =  2, 10, 18,  the  probability  of 
coexistence  of  a  secondary  link  is  plotted  in  Fig.  5  as  a 
function  of  2th-  While  the  Table  III  algorithm  behaves  - 
arguably-  satisfactorily,  it  can  be  observed  that  the  active 
interference  avoidance  algorithm  of  Table  IV  can  significantly 
enhance  the  opportunity  of  coexistence  of  a  secondary  link. 

Next,  under  the  same  primary  system  simulation  environ¬ 
ment,  we  examine  the  problem  of  cognitive  code-division 
channelization  for  multiple  secondary  links.  The  receiver 
output  SINR  requirement  for  each  secondary  link  is  set  at 
7 j  =  lOdB,  j  =  1,2,....  The  average  number  of  activated 
secondary  links  by  the  algorithm  in  Table  V  is  plotted  in  Fig. 
6  as  a  function  of  the  interference  threshold  2th  for  varying 
primary  user  population  K .  In  Fig.  7,  we  plot  the  average 
number  of  active  secondary  links  versus  the  number  of  primary 
users  for  varying  interference  threshold  values  2th.  The  exper¬ 
iments  demonstrate  that  the  developed  multiple  secondary  link 
channelization  method  can  effectively  produce  multiple  active 
secondary  links  and  enhance  spectrum  efficiency. 
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Fig.  5.  Probability  of  coexistence  of  secondary  link  versus  interference 
threshold  ^th- 


Fig.  6.  Average  number  of  activated  secondary  links  by  the  algorithm  in 
Table  V  versus  interference  threshold  for  varying  primary  user  population 
K. 

VII.  Conclusions 

We  considered  the  general  problem  of  establishing  a  sec¬ 
ondary  code-division  link  alongside  a  primary  code-division 
multiple-access  system.  We  first  developed  a  novel  iterative 
least  square  (ILS)  based  primary-system  identification  algo¬ 
rithm  which  can  blindly  detect  the  code  channels  utilized  by 
primary  users.  Then,  we  proposed  two  alternative  schemes, 
one  of  low  (passive  scheme)  and  one  of  moderate  (active 
scheme)  computational  complexity  that  optimize  transmitting 
power  and  binary  code-channel  allocation  of  the  secondary 
link  without  causing  “harmful”  interference  to  the  primary 
users.  At  the  same  time,  the  signal-to-interference-plus-noise 
ratio  (SINR)  of  the  secondary  link  at  the  output  of  the 
maximum  SINR  linear  receiver  is  no  less  than  a  certain 
threshold  to  meet  quality  of  service  (QoS)  requirements  for 
the  secondary  link.  Finally,  we  extended  the  channelization 
problem  to  multiple  code-division  secondary  links. 

Simulation  results  demonstrated  that  the  proposed  blind 
identification  algorithm  can  efficiently  and  effectively  detect 
primary  users’  code  channels  and  the  proposed  code-division 
channelization  methods  can  successfully  allow  secondary  links 
to  opportunistically  share  the  spectrum  without  causing  harm¬ 


Fig.  7.  Average  number  of  secondary  links  versus  number  of  primary  users 
K. 

ful  interference  to  primary  users. 

Cognitive  code-division  radios  combine  in  principle  the 
bandwidth  efficiency  characteristics  of  cognitive  operation 
and  code-division  multiple  accessing  and  are  expected  to 
find  a  place  in  future  communication  systems.  To  that 
extend,  the  developments  presented  in  this  paper  constitute 
an  early  contribution  that  can  be  helpful  in  benchmarking 
future  efforts.  Since  the  developed  herein  coexistence-mode 
cognitive  networks  do  not  require  any  prior  knowledge  about 
the  primary  networks  (such  as  signatures,  transmission  energy, 
channel  state  information),  the  secondary  networks  can  be 
deployed  transparently  without  any  modification/upgrade  of 
the  existing  primary  CDMA  infrastructure. 

Appendix  A 

Derivation  of  Iterative  Least-Squares  Procedure 
in  Table  I 

The  LS  cost  function  in  (15)  can  be  rewritten  as 
J  =  || Y  —  VB||f.  =  tr{YYff} -tr{YBffVff} 

-trjvBY^j +tr{vBBffVff}  (35) 

where  tr{-}  denotes  the  trace  of  a  matrix. 

For  a  given  B,  the  LS  optimal  estimate  of  V  can  be  obtain 
by  differentiating  the  cost  function  J  with  respect  to  Vn  and 
setting  the  outcome  equal  to  the  zero  matrix, 

8  J 

—jj  =  -YBff  +  V{BBh)  =  0  V  =  YBh(BBh)~1. 

(36) 

Next,  pretend  that  V  is  known  and  relax  the  domain  of  the 
symbol  information  matrix  to  the  complex  space,  B  £  CKxN . 
The  LS  optimal  estimate  of  B  £  CKxN  can  be  calculated 
again  by  differentiation 
8  T 

-—j  =  -VhY  +  YhVB  =  0  =>  B  =  (VhY)~1VhY. 
o 

(37) 

Finally,  we  project  (quantize)  the  complex-valued  LS  estimate 
B  to  the  binary  domain 

B  =  sgn  {93c  [(VffV)_1  V*Y]  }  .  (38) 
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